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ABSTRACT 

Context. Significant progress has been made in the last years in the understanding of the jet formation mechanism through a combi- 
nation of numerical simulations and analytical MHD models for outflows characterized by the symmetry of self- similarity. Analytical 
radially self-similar models successfully describe disk-winds, but need several improvements. In a previous article we introduced 
models of truncated jets from disks, i.e. evolved in time numerical simulations based on a radially self-similar MHD solution, but 
including the eff'ects of a finite radius of the jet-emitting disk and thus the outflow. 

Aims. These models need now to be compared with available observational data. A direct comparison of the results of combined 
analytical theoretical models and numerical simulations with observations has not been performed as yet. This is our main goal. 
Methods. In order to compare our models with observed jet widths inferred from recent optical images taken with the Hubble Space 
Telescope (HST) and ground-based adaptive optics (AO) observations, we use a new set of tools to create emission maps in diff'erent 
forbidden lines, from which we determine the jet width as the full- width half-maximum of the emission. 

Results. It is shown that the untruncated analytical disk outflow solution considered here cannot fit the small jet widths inferred by 
observations of several jets. Furthermore, various truncated disk- wind models are examined, whose extracted jet widths range from 
higher to lower values compared to the observations. Thus, we can fit the observed range of jet widths by tuning our models. 
Conclusions. We conclude that truncation is necessary to reproduce the observed jet widths and our simulations limit the possible 
range of truncation radii. We infer that the truncation radius, which is the radius on the disk mid-plane where the jet-emitting disk 
switches to a standard disk, must be between around 0. 1 up to about 1 AU in the observed sample for the considered disk- wind solu- 
tion. One disk- wind simulation with an inner truncation radius at about 0. 1 1 AU also shows potential for reproducing the observations, 
but a parameter study is needed. 

Key words. MHD — methods: numerical — ISM: jets and outflows — Stars: pre-main sequence, formation 



1, Introduction consistently treating the disk- jet system semi-analytically (e.g. 

^5 r—, — 1 1 -. lFerreirairi997l: ICasse & Ferreirall20QQah . or, b y self-consistently 

g Astrophysical jets and disks (|Livio| 12009D seem to be inter- treating numerically th e disk-jet system fe.g. lZanniet al.ll2007[ 

related, notably in young stellar objects (YSOs), where jet sig- [Tzeferacos et al ll2QQ9l) 






. natures are well correlated wit h the infrared excess and accretion 

■ rate o f the circumstellar disk ( Cabrit etaPl 19901: iHartigan et alJ The original lBlandford & Paynd (Il982h model, however, has 
[ |2004|) . Disks provide the plasma which is outflowing in the serious limitations for a needed meaningful comparison of its 
t, jets, while jets in turn provide the disk with the needed angu- predictions with observations. First, singularities exist at the jet 
C^ lar momentum re moval so that ac cretion onto the protostellar axis, the outflow is not asymptotically super-fast, and most im- 
object takes place ('Hartmann'2009). On the theoretical front, the portantly, an intrinsic scale in the disk is lacking with the result 
most widely accepted description of this accretion- ejection phe- that the jet formally extends to radial infinity, to mention just 
nomenon (Ferreira 2007) is based on the interaction of a large a few. First, the singularity at the axis can be easily taken care 
scale magnetic field with an accretion disk around the central of by numerical simulations extending the analy tical solutions 
object. Then, plasma is channeled and magneto-centrifugally close to this symmetry axis (iGracia et al.l l2006l GVT06 here- 
accelerated along the open magnetic field lines threading the after). Next, the outflow speed at large distances may be tuned 
accretion disk, as first described in Blandford & Payne (1982). to cross the corresponding limiting characteristic, with the result 
Several works have extended this study either by semi-analytic that the terminal wind solution is causally disconnected from the 
models using radially self- similar solutions of the full mag- disk and hence perturbations downstream of the super-fast tran- 
netohydrodynamics (M HD) equations with the disk treated as sition (as modified by self- similarity) ca nnot aflTect the whole 
a boundary condition (iVlahakis & Tsinganosl Il998l) , by self- structure of the steady disk- wind outflow (Vlahak is et al.ll200Q[ 

VOO hereafter) , a state which has al so been shown to be struc- 
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vious paper (IStute et alJ I2Q08L paper I hereafter), wherein we 
presented numerical simulations of truncated flows whose ini- 
tial conditions are based on analytical self- similar models. 

In order to test our truncated models, we will now apply 
our simulations to observations. In recent years, many NIR and 
optical data have become available exploring th e morphology 
and kinematics ofthe jet launching region (e.g. Doug ados et al. I 
i200Q: Rav etani2007l: IPougados 2008, and references therein). 
Hubble Space Telescope (HST) and adaptive optics (AO) ob- 
servations give access to the innermost regio ns of the wind , 
where the acceleration and collimation oc curs (IRav et al.ll 19961: 
iDougadosetaTlbOOOl: IWoitas et al.|[2002t iHartigan et al.l l2004). 
Because YSO jets emit in a number of ato mic (and molecular) 
lines, we used a set of tools described in iGracia et aP (|201Q|) 
to create emission maps in difl'erent forbidden lines which were 
used by other authors to extract the jet width from images. 
The observed jet widths will be compared with those extracted 
from our synthetic i mages. A similar st udy h as been done by 
ICabrit et al.1 (Il999l) . iGarcia et al.1 (l200ll) and iDougados et"aLl 
(200?) using a difl'erent se t of se mi-analytical self-similar disk- 
jet solutions from Ferreira (1 19971) and|Casse & Ferreira (2000a). 
Observed jet widths could be reproduced by manually truncating 
the solutions inside 0.07 AU and outside 1 AU, but the modifi- 
cation of flow streamlines induced by truncation was ignored. 

The remainder of the paper is organized as follows: we 
briefly review the initial setup of the numerical simulations in 
Sect. [2] and describe our procedure for the comparison with ob- 
servations in Sect. [S] The results of our studies are presented in 
Sects, m - [6j Finally, we conclude with the implications of the 
results in terms of the structure of the disk and the respective 
launching radii of the jets in YSOs. 



2. Initial model setup and numerical simulations 

This work is based on the results of our numerical simulations 
discussed in our paper I and two new mo dels. We solved the 
MHD equations with the PLUTO cod^H (Mign oneet"aDl2007h 
starting from an initial condition set according to a steady, ra- 
dially self-similar solution as described in VOO, hereafter la- 
belled ADO (analytical disk outflow solution, as in MOB), which 
crosses all three critical surfaces. At the symmetry axis, the ana- 
lytical solution was modified as described in GVT06 and M08. 
To study the influence of the truncation of the analytical so- 
lution, we divided our computational domain into a jet region 
and an external region, separated by a truncation field line Oftmnc- 
For lower values of the normalized magnetic flux function, i.e. 
(^ < ^tmnc - or conversely smaller cylindrical radii - our ini- 
tial conditions are fully determined by the solution of VOO and 
the modification of GVT06 and M08 close to the axis. In the 
outer region we modified all quantities and initialized them with 
another analytical solution, but with modified parameters. From 
VOO, one can show that if we start with an arbitrary MHD so- 
lution with the variables p, p, v, B, one can easily construct a 
second solution by using two free parameter^ Ai and A2, with 



V. Thus some 



p' = AiAIp,p' = AIp,B' =^2^andv' = A~^^^ 

or all quantities are scaled down in the external region depending 

on our choice of parameters. In Table [T] we give the parameters 



^ http://plutocode.to.astro.it/ 

^ Strictly speaking this is only true when all lengths are also scaled, 
but because only the gravity term explicitly depends on the length scale 
and in our case the gravitational force is small compared to the other 
forces, this slight inconsistency is unimportant (see paper I). 



of the models (some of which we studied in paper I as well as 
new ones) used in this study. 

For further details, we refer the reader to paper I. 

3. Comparison of synthetic emission runs with 
observations 

Numerical simulations and observations cannot be directly com- 
pared. While the former describe the plasma in terms of phys- 
ical quantities like density, pressure, magnetic field and veloc- 
ities, the latter observes only photon flux as a function of fre- 
quency. This comparison can be facilitated by means of con- 
structing synthetic maps. However, translating numerical simu- 
lations into such synthetic maps is very complicated. In general, 
the local emissivity is a function of temperature, electron density 
and density of the respective ion, as e.g., ^oii for singly ionized 
oxygen. The emissivities are integrated along a given line-of- 
sight and projected onto the plane of the sky producing an ideal 
synthetic emission map. Finally, real detectors distort this ideal 
map through their detector response, which needs to be taken 
into account. 

In this paper, we use the se t of tools OpenSESAMqj vO.l de- 
scribed in lGraciaet al.l ("2010") to produce synthetic observations 
from our simulations in difl'erent consecutive stages. The first 
stage approximates the ionization state of each atom by locally 
solving a chemical network under the assumption of local equi- 
librium. The second stage calculates the statistical equilibrium of 
level populations for each ion of interest as a function of temper- 
ature and density and yields the emissivity for individual tran- 
sitions of interest. Further stages take care of integration along 
the line-of-sight and projection. Finally, the ideal maps are con- 
volved with a Gaussian point- spread-function (PSF) to mimic 
the finite spatial resolution of a given instrument. These syn- 
thetic emission maps are then quantitatively analyzed with simi- 
lar techniques used on real observed maps. We refer to "runs" as 
runs of OpenSES AMe with diff'erent sets of scalings (see below) 
and to "models" as diff'erent MHD simulations with PLUTO 
(TablelB. 

We comp are the width of jets nieasured from HST and AO 
obser vations (iDougados et al. iboOOl: iRav et al.ll2007l: iDougadosI 
'2008', and references therein) with the width from synthetic 
emission maps calculated from our MHD models. We convolved 
the maps with a Gaussian with a full-width half-maximum 
(FWHM) of 15 AU (cr = 6.37 AU) throughout this paper, equiv- 
alent to HST's resolution of 0.1" at a distance of 150 pc. We use 
a sample of eight jets: DG Tau, HN Tau, CW Tau, UZ Tau E, 
RW Aur, HH34, HH30 and HL Tau (Fig. [B. In order to deter- 
mine the width of the jets in our models, we use a method which 
is as close as possible to that applied by the observers. We cre- 
ate convolved synthetic maps for the emission in the [SII] /16731 
and [01] /16300 lines for each numerical model and each run 
of OpenSES AMe and determine the jet width from the map's 
FWHM as a function of distance along the axis. 

3.1. Normalizations 

Throughout the first paper, we used only the dimensionless 
quantities in which PLUTO performs its calculations. In order 
to compare our results with observations, however, i.e. in order 
to run OpenSESAMe correctly, we have to scale them to phys- 
ical units by providing scaling factors for density po, pressure 
po, velocity vq, magnetic field strength ^o, a length scale Rq and 



http ://homepages . dias .ie/~j gracia/OpenSES AMe/ 
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Table 1. List of numerical science models 



Name 



Domain [Rq x Rq] Resolution Description 



= 0.4, external analytical solution Ai = \0^, A2 = 10""^ 
= 0.2, external analytical solution Ai = \(fi, A2 = 10~^ 
= 0.1, external analytical solution Ai = 10^, /I2 = 10~^ 
= 0.01, external analytical solution Ai = 10^, A2 = 10~^ 
= 0.001, external analytical solution Ai = 10^, A2 = 10~ 
= 0.4, external analytical solution Ai = 100, A2 = 0.1 
same as model SC2, but solutions are swapped 
Q^trunc = 0-4, external analytical solution /li = 1, /I2 = 0.1 
same as model SC4, but solutions are swapped 



^trunc [^0] 



model 
model 
model 
model 
model 
model 
model 
model 
model 



SCla 

SClb 

SClc 

SCld 

SCle 

SC2 

SC3 

SC4 

SC5 



[0,50] X 
[0,50] X 
[0,50] X 
[0,50] X 
[0,50] X 
[0,50] X 
[0,50] X 
[0,50] X 
[0,50] X 



[6,100] 
[6,100] 
[6,100] 
[6,100] 
[6,100] 
[6,100] 
[6,100] 
[6,100] 
[6,100] 



200 X 400 
200 X 400 
200 X 400 
200 X 400 
200 X 400 
200 X 400 
200 X 400 
200 X 400 
200 X 400 



<3^trunc 
<3^trunc 
Q^trunc 
<3^trunc 
<3^trunc 
Q^trunc 



model self 
model SClg 



[0,50] X [6,100] 
[0,10] X [6,20] 



200 X 400 
200 X 400 



Q^trunc = 0.0005, external analytical solution Ai = 10-^, /I2 = 10 
Q^trunc = 1 X 10~^, external analytical solution Ai = 10^, A2 = 10~^ 



5.375 
5.125 
4.875 
3.625 
2.625 
5.375 
5.375 
5.375 
5.375 



2.375 
0.575 



200 



^ 150 



100 



50 



DGTau o CWTau □ RW Aur o HL Tau --■ 
HN Tau + UZ Tau E HH 30 — HH 34 




^Sfooooooooo' 



100 200 300 400 500 600 
distance from jet source [AU] 



700 



800 



Fig. 1. Variation of jet width (FWHIM) derived from [SII] and 
[01] images as a function of distance from the source. Data 
points are from CFHT/PUEO and HST/STIS observations of DG 
Tau (diamonds), HN Tau (plus signs), CW Tau (squares), UZ 
Tau E (crosses), RW Aur (circles), HH 34 (one triangle), HH 30 
(blac k solid line) and HL Tau (red dashed line); data are taken 
from lRay et alJ (l2QQ7l) and references therein. 



a mass scale M. However, in terms of the normalizations used 
in the PLUTO code, only three of those quantities are indepen- 
dent. A possible choice is the mass of the central object, velocity 
scale and density scale, while the remaining factors are calcu- 
lated from these. 

Here we use three different "coordinate systems": i) the com- 
putational grid of cells with indices (ij) from (0,0) to (199,399), 
ii) the PLUTO domain from (0,6) to (50,100) and iii) the phys- 
ical scale of the jet in AU, which is simply the PLUTO domain 
multiplied with a length scale Rq. In the solution of VOO, the 
length scale Rq is connected to the mass of the central object and 
the velocity normalization via 



4 v^ Vkms V 



-2 



M 



0.5 M. 



(1) 



From the velocity and density normalization directly follow 
the normalizations for the magnetic field and pressure as 



Po = PoVq = 10 ^^gcm ^ s 



Po 



gem" 



M^l^'^ 



''^^^^'-''■'"^{w^fi^)- 



(3) 



The mass of the central object aff'ects only the length scale. 
The pressure and temperature of the jet and thus the synthetic 



emission maps are aff'ected only by the unit density and unit ve- 
locity. 

As typical jet velocities in YSOs we assumed values of 100, 
300, 600 and 1000 k m s~\ as typical ma sses of T Tauri stars 
0.2, 0.5 and 0.8 M© (Hartigan et al. "1995), and as jet number 
densities values of 125, 500, 1000 and 5 x 10"^ cm'^ We adopt 
the nomenclature for our runs as e.g. (^jet, Vjet, M) with ^jet in 
cm"^, Vjet in km s"^ and M in M©. 

In order to find the normalizations listed above, we needed to 
have typical values of jet density and velocity at a certain posi- 
tion (R, z) for a given mass M, and iteratively solved the Eq. ([T]). 
Details of this algorithm and also a graphical picture of this ap- 
proach are given in Appendix |A1 The normalizations vary from 
numerical model to numerical model, therefore we used the cor- 
responding diff'erent values for each model. 

As another constraint, we require that Rq is small enough that 
the FWHM of the Gaussian of 15 AU is sampled by a reason- 
able number of pixels. Because the resolution of our numerical 
simulations was four pixels per Rq, requiring at least two pixels 
per FWHM requires Rq < 30 AU. This limit highly reduces the 
number of runs, from 576 to 180. The only valid runs are 

- (• • -,600,0.2), (• • -,600,0.5), (• • -,1000,0.5), (• • -,1000,0.8) for 
models ADO, SCla-b, SC2, SC4 

-(•••,600,0.2), (•••,600,0.5), (•••,1000,0.2), (•••,1000,0.5), 
(• • •,1000,0.8) for models SClc-f 

- (•••,600,0.2), (•••,1000,0.2), (•••,1000,0.5), (•••,1000,0.8) 
for model SClg 

- (•• •,100,0.2) for model SC3 

- no valid runs for model SC5. 

The corresponding values of Rq are listed in Tabled 



3.2. The artefact of limb-brightening 

After plotting transverse intensity cuts through the convolved 
emission maps, it can be seen that the maximum emission does 
not come from the axis, but from a small shell at a finite radius 
(Fig. O left). This eff'ect of limb-brightening, which has been 
never observed up to now in real proto stellar jets, is a direct con- 
sequence of the density and pressure/temperature profiles in the 
analytical solution and is present in all our models and runs, the 
untruncated model ADO as well as our truncated models. Details 
of the physical reasons for this behavior are given in Appendix 

m 

Because the analytical solution is not well defined close to 
the axis and we had to interpolate it in our numerical simula- 
tions, and bec ause it is like ly that a stellar wind resides inside 
the disk wind (iMatsakos et al. 2009), we applied the following 
corrections before running OpenSESAMe: 
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Fig. 2. Transverse intensity cuts of the convolved [01] /16300 
emission at z = 100 AU (black, solid), normalized to its value 
on the axis, for model ADO and run (500,600,0.5); without cor- 
rections of density and temperature they show the effect of limb- 
brightening (top) which vanishes after the corrections (bottom); 
also plotted are cuts of the normalized hydrogen number density 
(long-dashed, green), of noil no (blue, dotted) and of the tem- 
perature (dashed, red, y-axis on the right-hand side). 7?o = 7.60 
AU. 

- we limited the temperature to 10^ K in the whole domain 

- we limited the density around the axis by setting the density 
inside 1 FWHM to its value at 1 FWHM for each z. 

This value of T is typical o f that deduced from anal- 
yses of line ratios in jets (e.g. iBacciotti & EisloffeJI 119991: 
iLavalley-Fouquet et al.ll200d) . The artefact of limb-brightening 



is removed by limiting T . For numerical reasons, the density 
shows a steep increase close to the jet axis. This artefact has 
to be corrected to obtain jetwidths above 1 FWHM. 

4. Do we have to truncate the disk? 

Before we present the implications of truncating the analytical 
solutio n below, we have to compare the unmodified ADO so- 
lution (IVlahakis et al.ll200U ) with observations to show that the 
truncation is indeed needed. Figure[3] shows the typical output of 
OpenSESAMe: cuts of the density and temperature, the electron 
and SII and 01 ion densities and synthetic emission maps of the 
[SII] ^6731 and [01] ^6300 lines, all for model ADO and run 
(500,600,0.5). The electron density shows a steep cut-off' along 
a field line anchored at about 5 7?o where the temperature drops 
below 10"^ K, while the [SII] and [01] emission maps show a 
somewhat shallower cut-off', which is located at a larger radius 
(about 8 7?o). Synthetic emission maps of the [01] /16300 line for 
all models and runs are shown in Appendix O 

The extracted jet widths in normalized 7?o units are plotted in 
Fig. m In each plot we combine runs with equal velocities and 
masses, i.e. only varying the density, which result in identical jet 
widths as expected. Because we extracted the jet width from a 
ratio of intensities by using the FWHM (we divided the maps 
by the intensity on the axis and check where the ratio is 0.5), the 
factor p^ cancels out. 

Much larger changes are present when we compare runs 
where po is identical, but vq varies (top right and bottom left 
plot) and when we look for the inffuence of the mass (plots in 
each row). The jet widths derived from synthetic [01] are iden- 
tical to those derived from synthetic [SII] images, thus we show 
only results based on [01] images. 

The jet widths in Rq increase with increasing velocity. A 
higher mass reduces the normalized jet width considerably. In 
run (500,600,0.5), as seen in Fig.O the extracted jet width does 
not follow any directly apparent feature in the emissivity maps 
nor a specific contour line. 

After rescaling the jet widths with the appropriate values of 
7?o, we compared the modeled and observed jet widths in AU 
(Fig.O. The jet widths are now rescaled with 7?o, which is pro- 
portional to the mass and anti-proportional to vq (eq.[T]). Because 
/(/) = Vjet/vo is almost constant in our range of /?o, Vjet and vq 
can be interchanged. This behavior of 7?o is now dominant, thus 
the jet width increases with increasing mass and decreases with 
increasing 7^r velocity, although the second eff'ect seems to be of 
minor importance. 

The untruncated ADO model on average gives too large jet 
widths compared to the observations of T Tauri jets. The first 
five data points within 120 AU from the jet source in DG Tau 
can be best approximated with run (500,600,0.2), unfortunately 
our model does not provide results farther out for such a small 
Rq. At distances above 200 AU, the run (500,600,0.5) gives a jet 
width in the observed range but only after the first bump, w hich 
is intrinsic for our analytical model. iHartigan et al.l (1 19951) give 
a mass of DG Tau of 0.67 M©, i.e. higher than in both models. 
Although this diff'erence in mass might not be meaningful, we 
conclude that to be able to reproduce all jets in our sample, we 
need an additional eff'ect which reduces the derived jet width. 
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Fig. 3. Typical output of OpenSESAMe: cuts of the density and temperature, the electron and SII and 01 ion densities and synthetic 
emission maps of the [SII] ^673 1 and [01] /16300 lines for model ADO and run (500,600,0.5). Rq = 7.60 AU. 
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Fig. 4. Dimensionless jet width derived from synthetic [01] images as a function of distance from the source in model ADO. 
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Fig. 5. Jet widths in AU derived from synthetic [01] images as 
a function of distance from the source in model ADO; overlaid 
are the data points of Fig. \T\ 



5. Effects of outer truncation 

First we investigate the effects of outer truncation. Convolved 
synthetic maps for the emission in [01] /16300 for some numeri- 
cal models and run (500,1000,0.5) are given in Fig.[6j Truncation 
leads to collimation of the emission region with respect to the 
model ADO without any truncation. 

Again we extracted the jet width from emission maps like 
these. The resulting widths derived from the synthetic [01] im- 
ages and scaled to AU are presented in Fig. [71 We found sim- 
ilarities in behavior in the truncated models to the untruncated 
model ADO. The jet widths show again no dependency on 
the density, as described for model ADO in the previous sec- 
tion. Surprisingly, in models SCla-c, SC2 and SC4 the runs 
(500,600,0.2) and (500,1000,0.5) and also (500, 1000,0.8) lead to 
almost similar physical jet widths. The first two also almost coin- 
cide in models SCld-e. As in model ADO, also in the truncated 
models the run (500,600,0.5) has the smallest jet widths (after 
the first bump). In principle, we can reproduce even smaller val- 
ues than the observed ones. 



6. Effects of inner truncation 

In paper I we also performed numerical simulations, in which 
we truncated the analytical solution in the interior, i.e. at an 
inner truncation radius. The physical picture behind this sce- 
nario is a stellar magnetosphere truncating the jet-emitting disk. 
We showed that inner truncation leads to a decrease of the 
jet radius and compression of the material in the inner region. 
Unfortunately, only one run met our scaling requirements (Sect. 
13.1b : model SC3 and run (• • •, 100, 0.2). For this model, a con- 
volved synthetic map for the emission in [01] /16300 is given in 
Fig.m 

After rescaling the derived jet width to AU, we found an al- 
most constant width in the range of the observed values. Note 
that our model does not provide results farther out than 100 AU 
due to a small 7?o. 
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Fig. 8. Synthetic emission maps of the [01] /16300 line, con- 
volved with a Gaussian PSF with a FWHM of 15 AU, for model 
SC3 and run (• • •, 100, 0.2). Rq = 0.97 AU. 
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Fig. 9. Jet widths in AU derived from synthetic [01] images as a 
function of distance from the source in model SC3 and for run 
(500,100,0.2); overlaid are the data points of Fig. 1. 



We found that the untruncated model ADO of lVlahakis et al.l 
(2000) cannot account for the small jet widths found in recent 
optical images taken with HST and AO. The density normaliza- 
tion is not important for the resulting measured jet width as long 
as we are far from the critical regime. 

We investigated different effects for reducing the deriving jet 
width: by imposing an outer radius of the launching region of 
the underlying accreting disk and thus also of the outflow on the 
observable structure of the jet and by imposing an inner radius of 
the underlying accretion disk due to interactions with the stellar 
magnetosphere. 



7. Discussion 

We studied the jet widths derived from synthetic emission maps 
in different forbidden lines as the full-width half-maximum of 
the emission. 



7.1. Outer truncation 

We created synthe tic images based on our simulations of trun- 
cated disk winds ( Stute et al.ll2008l) as well as new simulations 
and found that the extracted jet widths in the truncated models 
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Fig. 6. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for run 
(500,1000,0.5). Top row: models ADO (Rq = 1.04 AU), SCla (Rq = 1.98 AU), SClb (Rq = 1.98 AU), SClc (Rq = 2.19 AU); 
bottom row: models SCld (Rq = 2.61 AU), SCle (Rq = 4.22 AU), SClf (Rq = 4.86 AU), SClg (Rq = 7.23 AU). 



decrease for models SCla-lg, compared to those of the untrun- 
cated model ADO, as naively expected. 

In the present paradigm, jets are emitted only by the inner 
part of the disk. Hence in the other parts the disk can be de- 
scribed by a standard accretion disk (SAD), in th e inner parts 
by a jet-emitting disk (JED). Anderson et alj (l2003h showed that 
one can estimate the launching region as 






\2/3 



UoOkms-V \0. 



lOkms-i^ 

1/3 



M 



5M« 



This transition was constrained observationally with mea- 
sured jet rotation velocities and using the equation above and 
radii of the o rder of 0.1-1 AU are used in several theoretical 
studies as e.g. lCombet & Ferre ira (2008). 

Our results can be used to infer the "real" value of the trun- 
cation radius 7?tmnc in the observed sample of jets and interpret 
it as the transition radius of the JED to the SAD, assuming the 
specific model of VOO applies. 

At the lower boundary in our simulations, the truncation radii 
are given in Table[TJ They vary from 5.375 Rq in model SCla to 
0.575 Rq in model SClg. However, these radii are set at z = 67?o 
(the lower boundary), not in the equatorial plane. Those can be 
calculated by extrapolating the field line, i.e. 



R^ 



tmnc lz=0 



= /?, 



tmnc \z-6 



Gjnll) 



(4) 



with ^tmnc = arctan(7?tmnclz=6/6) and G taken from the analytical 
solution of VOO. This gives the following results: 

- in models SCla, SC2-4 7?tmnclz=o = 0.1 1316 7?o, 

- in model SClb 7?tmnclz=o = 0.09099 7?o, 

- in model SClc 7?tmnclz=o = 0.07275 /?o, 



- in model SCld /?tmnclz=o = 0.0213 7?o, 

- in model SCle 7?tmnclz=o = 0.00722 7?o, 

- in model SClf 7?tmnclz=o = 0.00546 7?o and 

- in model SClg 7?tmnclz=o = 0.00069 7?o. 

The last step is a multiplication with 7?o of the best-fit model 
(TableO. The resulting values of 7?tmnclz=o in AU for each model 
and run are given in Tabled 

We found best-fit models for the jets in the observed sample. 
Note that we always ignored the first bump in the synthetic jet 
widths and focussed on larger distances from the source: 

- The observed mass of DG Tau (diamonds in Fig. [7]) is 0.67 
Mq (Hartigan et al.l 1 1995b . therefore we have to focus on 
the runs (500,1000,0.8), and perhaps also runs (500,600,0.5) 
and (500,1000,0.5). The best-fit model is between ADO and 
SCla, thus the truncation radius is larger than 0.22 AU. 

- HN T au (plus signs) has a mass of 0.72 M© (iHartigan et al.l 
119951) . thus again the runs (500,1000,0.8) are favored. 
Because we ignored the first bump, we interpolated the jet 
shape at larger distances and found a best-fit model between 
ADO and SCla, thus the truncation radius is again larger 
than 0.34 AU. However, this result is highly uncertain. 

- The mass of CW Tau (squares) is th e highest in our 
sample, 1.03 M© (IHartigan et al.l 1 19951) . Using the runs 
(500,1000,0.8), the best-fit model is SClb or SClc. The trun- 
cation radius is thus between 0.25 - 0.3 AU. 

- UZ Tau E (crosse s) has the lowest m ass in our sam- 
ple, only 0.18 M© ("Hartigan et al."l995'), we use the runs 
(500,600,0.2). Again we had to interpolate the jet width from 
larger distances and choose model SCla as best-fit model. 
The truncation radius is about 0.26 AU. 

- The measured mass of RW Aur (circles) is 0.85 M© 
([Hartigan et al.J Il995l) . thus we have to focus on runs 
(500,1000,0.8). We need a very high degree of truncation 
as in models SCle-g, thus a truncation radius of the order of 
0.04 AU. 
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Table 2. Rq in AU of the different OpenSES AMe runs 



model 


(• • -,600,0.2) 


(• • ^,600,0.5) 


(•••,1000,0.2) ( 


•••,1000,0.5) 


(•••,1000,0.8) 


ADO 


1.20 


7.60 


— 


1.04 


2.47 


SCla 


2.32 


7.60 


— 


1.98 


2.99 


SClb 


2.17 


7.59 


— 


1.98 


3.33 


SClc 


2.40 


7.60 


0.98 


2.19 


3.47 


SCld 


2.91 


8.87 


1.34 


2.61 


4.10 


SCle 


4.54 


13.96 


2.06 


4.22 


5.86 


self 


5.29 


14.40 


2.17 


4.86 


6.91 


SClg 


8.03 


— 


3.04 


7.23 


13.69 


SC2 


2.31 


7.60 


— 


1.96 


2.98 


SC4 


2.05 


7.60 


— 


2.08 


2.97 


model 


(•••,100,0.2) 


SC3 


0.97 



Table 3. Truncation radius at the equator /?tmnclz=o in AU of the different OpenSES AMe runs 



model 


(• • ^,600,0.2) 


(• • ^,600,0.5) 


(•••,1000,0.2) 


^••,1000,0.5) 


(•••,1000,0.8) 


SCla 


0.263 


0.860 


— 


0.224 


0.338 


SClb 


0.197 


0.691 


— 


0.180 


0.303 


SClc 


0.175 


0.553 


0.071 


0.159 


0.252 


SCld 


0.062 


0.189 


0.029 


0.056 


0.087 


SCle 


0.033 


0.101 


0.015 


0.030 


0.042 


self 


0.029 


0.079 


0.012 


0.027 


0.038 


SClg 


0.006 


— 


0.002 


0.005 


0.009 


SC2 


0.261 


0.860 


— 


0.222 


0.337 


SC4 


0.232 


0.860 


— 


0.235 


0.336 


model 


(•••,100,0.2) 


SC3 


0.11 



7.2. Inner truncation 

Because the jet-emitting accretion disk is thought to be truncated 
by the stellar magneto sphere, we also investigated our models 
in which the analytical solution is truncated at an inner radius. 
We found that inner truncation can also reduce the extracted jet 
widths. The jet width from our model SC3 and run (• • -,100,0.2) 
is about 30 AU for the inner 100 AU of the jet. This is well in 
the observed range of 15-45 AU. Because we did not vary the 
inner truncation radius, we can only claim - based on our results 
of paper I - that the larger the truncation radius, the higher is the 
compression of the resulting jet. In our models we chose a radius 
of 5.375 7?o at z = 6. For model SC3 and run (- - -,100,0.2), this 
corresponds to an inner truncation radius of 0.11 AU, which is 
only slig htly higher th an the inner hole in T Tauri disks (0.02- 
0.07 AU. lNajita et al. ||2Q0 7). A parameter study varying the in- 
ner truncation radius is needed. 



8. Conclusions 

We showed as a proof of concept that jet widths derived from 
numerical simulations extending analytical MHD jet formation 
models can be very helpful for understanding recently observed 
jet widths from observations with adaptive optics and space tele- 
scopes. However, further aspects have to be investigated in more 
detail. 

An intrinsic feature in all our models with outer truncation is 
the first bump in the extracted jet width, which complicates the 
comparison of synthetic and observed jet widths. Only for DG 
Tau and CW Tau, we could unambiguously find models with 
outer truncation which fitted the observed jet widths at larger 
distances. For HN Tau and UZ Tau E, we have no observed 
jet widths at larger distances, only in the region which is con- 



taminated by the bump. RW Aur shows very small observed jet 
widths at all scales, which cannot be reproduced by any of our 
models and runs. The derived truncation radii (> 0.25 AU) are 
several times larger than the inner radius of the gaseous disk in 
T Tauri stars (0.02-0.07 AU, Najita et al. 2007), thus our use of 
a self- similar disk- wind solution is consistent. 

The model with inner truncation gives a synthetic jet width, 
which is constant and is within the range of observed widths. 
Another advantage of this model is that its flow velocities of the 
order of 100 km s"^ are closer to observed values. All models 
with outer truncation have very high velocities (> 600 km s"^ at 
our scaling point of 7?jet = 10 AU and Zjet = 100 AU), several 
times higher than in that with inner truncation. 

Because we originally focused only on the effect of outer 
truncation, we kept the inner truncation radius constant in our 
models SC3 and SC5. In order to further explore the ability of 
inner truncation to reproduce the observations, we have to vary 
the inner truncation radius in a parameter study. Naively, one 
would expect that the derived jet widths are not decreasing for 
decreasing truncation radii as for the outer truncation, but are 
increasing. This, however, has to be tested with new simulations. 

In our study, we assumed an inclination of 90° of the jet, 
thus projection effects may slightly change our results. We will 
pre sent the results of su ch a study in a forthcoming paper. 

'Garci a et aP (l200ll) showed in their Fig. 1 that the measured 
jet width s are mainly ch aracterized by the ejection index ^, de- 
fined by Ferreiral (1 19971) . This is related to the model parameter 
X of the solution of VOO, ^ = 2 (x - 3/4), and because in our sim- 
ulations X = 0.75, we get f = , which is intrinsic for a standard 
disk. In the solution of' Ferreiral (11997 V ^ also controls the open- 
ing of the field lines , because it i s conn ected to the lever arm A 
by A = 1 + 1 /(2 ^). [G^cia et al.l (l200ll) favored cold solutions 
with values of A between 50-70. If heating along streamlines is 
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allowed, the relation is b roken and also warm solutions of e.g. 
ICasse & Ferreiral (l2QQQbh with smaller A values of 8 and open- 
ing of streamlines (the ratio of the maximum radius to the initial 
launch radius) of about 30 can reproduce the observation. In the 
solution of VOO, we have A = G(7i/2)~^ ^ 39 and the maxi- 
mum opening Gjnax/G(7i/2) ^893. The use of other solutions 
will therefore highly influence our results in terms of the outer 
truncation radius. 



Appendix A: Details of the normalization algorithm 

In order to find the reference values for the distance Rq, pres- 
sure po and magnetic field ^o Hsted in Eqs. O-©, we require 
some typical values of jet density and velocity at a certain po- 
sition - we define the jet as the region at the physical position 
(^jet, ^jet) = (10, 100) AU - and solve iteratively Eq. ©: 

1 . we start the iteration with the choice Ro = I AU, namely in 
the grid cefl (40,398) 

2. we find vq by dividing the required jet velocity by the veloc- 
ity in PLUTO units inside the grid cell 

3. we calculate Rq with equation ([T]) 

4. we calculate the physical position of the grid cell using this 
new/?o 

5. if R (z) of the grid cefl is larger than 10 AU (100 AU), we 
move to another grid cell 10 % closer to the jet axis (equa- 
torial plane); if R (z) of the grid cell is smaller than 10 AU 
(100 AU), we move to another grid cell 10 % further away 

6. we stop the iteration when the grid cell does not change any- 
more between two steps. 

The same results, perhaps more intuitively, can be achieved 
using the following graphical picture. We assume that only grid 
cells along a cone with indices (/, 10 /) are important. Using Eq. 
([T]) and the relation / = 4R/Ro, we have to search for the inter- 
section point 



Vjet = 600 km s-\ M = 0.5 M© 



/(0 = v(/,10/)=^=4Vjet 

Vq 



R^ 



jet 1 . .. 

— F = S(i) • 



CM V/ 



(A.1) 



Figure lAJl shows /(/) for all models and g(i) for Vjet = 600 km 
s~^ (vjet = 100 km s"^ for comparison with models SC3 and 
SC5) and M = 0.5 M©. After an initial drop close to the origin, 
/(/) converges to an almost constant normalized velocity for / > 
7 (i.e. Ro < 5.7 AU for 7?jet = 10 AU) in the models ADO, SCla- 
c and SC5 and for / > 20 (i.e. Ro <2 AU) in the others. Its value 
increases from the untruncated model ADO (100) monotonically 
with increasing degree of truncation (150 in SCla up to 300 in 
SClg). The models SC2 and SC4 have values close to that of 
model SCla. 

Reducing the required mass or increasing the required jet 
velocity shifts the function g(i) upwards and thus increases / of 
the intersection point and decreases Ro, until g(i) does not in- 
tersect with /(/) anymore inside the computational domain. For 
example, we could not find acceptable sets of normalizations for 
the runs (• • -,1000,0.2) in models ADO, SCla-b, SC2 and SC4; 
in the models SClc-g, normalizations were found for all runs. 
However, most of these runs were later excluded due to our re- 
quirement ofRo < 30 AU. 

In the inner truncation models SC3 and SC5, the final value 
of /(/) is much lower, only about 15, therefore both functions 
only intersect for small velocities. Again, most of the runs found 
were later excluded due to our requirement of Ro < 30 AU. 
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Fig. A.l. Shape of the functions /(/) for all models and of the 
function g(i) (thick curve) for Vjet = 600 km s"^ (100 km s"^ 
for the models with inner truncation in the bottom panel) and 

M = 0.5 Mo. 



Appendix B: Reasons for the artefact of 
limb-brightening 

If the number density is below the critical density, the emissivity 
of e.g. the [SII] /16731 line can be written as 



€(T) = ^e ^SII f(T) = ^e ^s g(T) f(T) , 



(B.l) 



with Hq and ^sii the electron and ion density, respectively. The 
latter is a function of the total sulfur density ^s and the ionization 
fraction g(T) = ^sii/^s- f(T) usually has the form 



f(T) 



T'^ exp - 



D/T 



(B.2) 



with a powerlaw and an exponential cut-off'. One can show that 
6(T) has a maximum at temperatures of about 4 x 10^ K. For 
lower temperatures, f(T) decreases steeply. For temperatures 
higher than about 6 x 10^ K, g(T) decreases steeply. The latter 
cut-off' is due to the assumption of ionization equilibrium used in 
OpenSES AMe for calculating the level population and chemical 
composition, while the former is due to the exponential temper- 
ature dependence in Eq. (IB.2I) , intrinsic to collisional excitation 
of a level of upper energy D. In eq. dB.lb , Hq is also a very steep 
function of T in ionization equilibrium. 

Thus the main emission comes only from a thin shell of the 
jet, in which the temperature is exactly in the range of (4-6)x 10^ 
K. By changing the temperature normalization, we can move the 
emitting shell on the axis, but the resulting jet widths are of the 
order of the FWHM of the Gaussian PSF. Therefore we have to 
have a constant temperature profile towards the axis. 



10 Matthias Stute et al.: Comparison of synthetic maps from truncated jet-formation models with YSO jet observations 

Appendix C: Emission maps of the [01] /16300 line 
for all models and runs 

Here we present the emission maps of the [01] /16300 line for all 
models and runs (500, •••,•••) for the sake of completeness. 
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Fig. 7. Jet widths in AU derived from synthetic [01] images as a function of distance from the source in models ADO and SCla- 
SC2 and SC4 and for runs (500,600,0.2), (500,600,0.5), (500,1000,0.5), (500,1000,0.8); overlaid are the data points of Fig. 1. 
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Fig. C.l. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model ADO 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.5), (500, 1000, 0.8). 
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Fig. C.2. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SCI a 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.5), (500, 1000, 0.8). 
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Fig. C.3. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SClb 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.5), (500, 1000, 0.8). 








Fig. C.4. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SClc 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.2), (500, 1000, 0.5), (500, 1000, 0.8). 
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Fig. C.5. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SCld 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.2), (500, 1000, 0.5), (500, 1000, 0.8). 
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Fig. C.6. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SCle 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.2), (500, 1000, 0.5), (500, 1000, 0.8). 








Fig. C.7. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SClf 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.2), (500, 1000, 0.5), (500, 1000, 0.8). 
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Fig. C.8. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SClg 
and runs (500, 600, 0.2), (500, 1000, 0.2), (500, 1000, 0.5), (500, 1000, 0.8). 
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Fig. C.9. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SC2 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.5), (500, 1000, 0.8). 
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Fig. C.IO. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SC3 
and run (500, 100,0.2). 
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Fig. C.ll. Synthetic emission maps of the [01] /16300 line, convolved with a Gaussian PSF with a FWHM of 15 AU, for model SC4 
and runs (500, 600, 0.2), (500, 600, 0.5), (500, 1000, 0.5), (500, 1000, 0.8). 



